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MULTISCALE TECHNIQUES FOR PARABOLIC EQUATIONS 


AXEL MALQVIST 1 ’ 2 AND ANNA PERSSON 1 


Abstract. We use the local orthogonal decomposition technique introduced 
in na to derive a generalized finite element method for linear and semilinear 
parabolic equations with spatial multiscale diffusion coefficient. We consider 
nonsmooth initial data and a backward Euler scheme for the temporal dis¬ 
cretization. Optimal order convergence rate, depending only on the contrast, 
but not on the variations in the diffusion coefficient, is proven in the Loo (U)- 
norm. We present numerical examples, which confirm our theoretical findings. 


1. Introduction 

In this paper we study numerical solutions to a parabolic equation with a highly 
varying diffusion coefficient. These equations appear, for instance, when modeling 
physical behavior in a composite material or a porous medium. Such problems are 
often referred to as multiscale problems. 

Convergence of optimal order of classical finite element methods (FEMs) based 
on continuous piecewise polynomials relies on at least spatial H 2 -regularity. More 
precisely, for piecewise linear polynomials, the error bound depends on ||u||# 2 , 
where ||u||jy 2 ~ e -1 if the diffusion coefficient varies on a scale of e. Thus, the 
mesh width h must fulfill h < e to achieve convergence. However, this is not 
computationally feasible in many applications. To overcome this issue, several 
numerical methods have been proposed, see, for example, 0, 0, i, na, na, 
na, and references therein. In particular, j! 6 | and na consider linear parabolic 
equations. 

In na a generalized finite element method (GFEM) was introduced and conver¬ 
gence of optimal order was proven for elliptic multiscale equations. The method 
builds on ideas from the variational multiscale method mm), which is based on 
a decomposition of the solution space into a (coarse) finite dimensional space and 
a residual space for the fine scales. The method in na. often referred to as local 
orthogonal decomposition, constructs a generalized finite element space where the 
basis functions contain information from the diffusion coefficient and have support 
on small vertex patches. With this approach, convergence of optimal order can be 
proved for an arbitrary positive and bounded diffusion coefficient. Restrictive as¬ 
sumptions such as periodicity of the coefficients or scale separation are not needed. 
Some recent works (0,0, m) show how this method can be applied to boundary 
value problems, eigenvalue problems, and semilinear elliptic equations. There has 
also been some recent work on the linear wave equation |Ij. 

In this paper we apply the technique introduced in Em to parabolic equations 
with multiscale diffusion coefficients. For the discretization of the temporal domain 
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we use the backward Euler scheme. Using tools from classical finite element theory 
for parabolic equations, see, e.g, ca, na, hhi, and references therein, we prove con¬ 
vergence of optimal order in the L 00 (L 2 )-norm for linear and semilinear equations 
under minimal regularity assumptions and nonsmooth initial data. The analysis is 
completed with numerical examples that support our theoretical findings. 

In Section [2] we describe the problem formulation and the assumptions needed 
to achieve sufficient regularity of the solution. Section [3] describes the numerical 
approximation and presents the resulting GFEM. In Section [4] we prove error es¬ 
timates and in Section [5] we extend the results to semilinear parabolic equations. 
Finally, in Section [6] we present some numerical examples. 


2. Problem formulation 


We consider the parabolic problem 

u - V • ( AVu ) = /, 
(2.1) u = 0, 

u(-,0) = u 0 , 


in n x (0, T], 
on dfl x (0, T], 
in f2, 


where T > 0 and U is a bounded polygonal/polyhedral domain in R d , d < 3. We 
assume A = A(x) and / = f(x,t), that is, the coefficient matrix A does not depend 
on the time variable. 

We let IJ 1 (U) denote the classical Sobolev space with norm 

IMI/rqn) = IMIz, 2 (fi) + ll^ w lli 2 (fi) 


and V = -ffo(U) the space of functions in I/ 1 (U) that vanishes on <9U. We use 
H~ 1 ( fl) = V* to denote the dual space to V. Furthermore, we use the notation 
L p (0,T; X) for the Bochner space with finite norm 

/•t - l /p 


F||l p (0,T;A-) 


IX 


dt 


1 < p < 00, 




(0,T;X) — ess sup 

0 <t<T 


\V X, 


where X is a Banach space equipped with norm || • ||_y- Here v £ H x (0,T; X) 
means v,v £ £ 2 ( 0 , T;X). The dependence on the interval [0,T] and the domain 
f2 is frequently suppressed and we write, for instance, £ 2 (^ 2 ) for ^ 2 ( 0 , T; L 2 (H)). 
Finally, we abbreviate the L 2 -norm || • || := || • ||l 2 (o) and the energy norm, |||'||| := 
||H 1/2v -||. 

To ensure existence, uniqueness, and sufficient regularity, we make the following 
assumptions on the data. 


Assumptions. We assume 

(Al) A £ L 00 ( 12 , K dxd ), symmetric, and 

n • r • r A(x)v • v 

0 < a := ess ml mi -, 

zen t>eR d \{o} v ■ v 

n A(x)v ■ v 

00 > p := ess sup sup - 

xen u£tt<*\{o} v ’ v 

(A2) uo £ L 2 , 

(A3) f,f£L 00 (L 2 ). 
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Throughout this work C denotes constants that may depend on the bounds a 


and /3 (often through the contrast /3/a), the shape regularity parameter 7 (3.1) of 


the mesh, the final time T, and the size of the domain f2, but not on the mesh size 
parameters nor the derivatives of the coefficients in A. The fact that the constant 
does not depend on the derivatives of A is crucial, since these (if they exist) are large 
for the problems of interest. This is sometimes also noted as C being independent 
of the variations of A. 

We now formulate the variational form of problem (2.1). Find u(-,t) £ V such 
that u(-, 0 ) = uq and 


( 2 . 2 ) 


(u, v ) + a(u, v ) = (/, v), Vv <E V, t £ (0, T\, 


where (u,v) = f^ uv and a(u,v ) = (AVu,Vv). 

The following theorem states existence and uniqueness for (2.2). The proof is 
based on Galerkin approximations, see, e.g., [ 6 j and [TO] . 

Theorem 2.1. A ssume that (Al), (A2), and (A3) holds. Then there exists a 
unique solution u to (2.2) such that u £ £ 2 ( 0 ,T;U) and u £ £ 2 ( 0 ,T; H -1 ). 


3. NuMERCIAL APPROXIMATION 


In this section we describe the local orthogonal decomposition presented in [14] 


to define a generalized finite element method for the multiscale problem ( 2 . 2 ). 

First we introduce some notation. Let { Th}h>o and {Tn}H>h be families of 
shape regular triangulations of fl where hx := diam(A’), for K £ 7 ~h, and Hx ■= 
diam(A'), for I\ £ Th ■ We also define H := m.&xxeT H Hk and h := maxx^T h h k ■ 
Furthermore, we let 7 > 0 denote the shape regularity parameter of the mesh Th ; 


(3.1) 


7 := max tr-. 
K&Th 


with jx '■= 


diam B x 
diarn K 


for K £ Th, 


where Bx is the largest ball contained in K. 

Now define the classical piecewise affine finite element spaces 


Vh = {z' £ C'(fl) : v = 0 on BQ,v\t is a polynomial of degree < l,VAf £ Th }, 
Vh = {'C £ C(£i) : v = 0 on d£t,v \t is a polynomial of degree < 1 , V/i £ Th}- 


We let A f denote the interior nodes of Vh and ip x the corresponding hat function for 
x £ A/", such that spand^}^^) = Vh ■ We further assume that Th is a refinement 
of Th, such that Vh Q Vh. Finally, we also need the finite element mesh Th of f 1 
to be of a form such that the L 2 -pi'ojection Ph onto the finite element space Vh is 
stable in H 1 - norm, see, e.g., [3], and the references therein. 

To discretize in time we introduce the uniform discretization 


(3.2) 0 = to < ti < ... < tx = T, where t n — t n - 1 = r. 

Let U n be the approximation of u(t) at time t = t n and denote /„ := f{t n ). Using 
the notation d t U n = (U n — C/„_i)/r we now formulate the classical backward Euler 
FEM; find U n £ Vh such that 

(3.3) (B t U n , v ) + a(U n , v ) = (/„, v), \/v £ V h , 

for n = 1,..., N and Uq £ Vh is some approximation of uq. For example, one could 
choose U 0 = P h u 0 , where Ph is the i^-projection onto 14. We also define the 
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operator A h : V h -> V h by 


(3.4) 


(AhV, w) = a(v, w), W, w G V h - 


The convergence of the classical finite element approximation (3.31 depends on 
||D 2 m||, where D 2 denotes the second order derivatives. If the diffusion coefficient 
A oscillates on a scale of e we have ||.D 2 w|| ~ e -1 . Indeed, defining A = —V • AX7, 
elliptic regularity gives 

||D 2 w|| < Ci||Au|| < C 2 \\AAu\\ < C 2 ||V • HVu - VH • Vu|| 

< C 2 (||^«|| + || VA ■ V«||) < C 2 (||>1«|| + Ca|| Vu||) < C 3 ( 1 + C a )\\Au\\, 

where C A is a constant that depends on the derivatives (variations) of A. This 
inequality is sharp in the sense that An and VH ■ V« does not cancel in general. 
The total error is thus || u(t n ) — U n \\ ~ (r + (h/e) 2 ), which is small only if h < e. 

The purpose of the method described in this paper is to find an approximate 
solution, let us denote it by U for now, in some space V C 14, such that dimt^ = 
dim Vh, for H > h, and the error \\U n — U n \\ < CH 2 . Here C is independent of 
the variations in A and U n is less expensive to compute than U n . The total error 
is then the sum of two terms 

I Ht n ) - U n \\ < || u(t n ) - U n \\ + || U n - U n ||, 

where the first term is the error due to the standard FEM approximation with 
backward Euler discretization in time. This is small if h is chosen small enough, 
that is, if h resolves the variations of A. Hence, we think of h > 0 as fix and 
appropriately chosen. Our aim is now to analyze the error \\U n — U n |. 

We emphasize that V = Vh is not sufficient. The total error would in this case 
be || u(t n ) — U n \\ ~ (r + (H/e) 2 ), which is small only if H < e. 

The next theorem states some regularity results for (3.3). 

Theorem 3.1. Assume that (Al), (A2), and (AS) holds. Then, for 1 < n < N, 
there exists a unique solution U n to (2.2) such that U n G 14. Furthermore, if 
U 0 = 0 , then we have the bound 

(3-5) \\dtU n \\<Cm L M + \\f\\ LM ), 

and, if f = 0, then 

(3.6) \\d t U n || < Ct-'WUoW, n > 1, ||^Cg| < Ctf 2 \\U 0 \\, n> 2, 

where C depends on a and T, but not on the variations of A. 


Proof. From (3.3) it follows for n > 2 that 

(d t d t U n , v) + a(8 t U n , v) = ( 8 t f n , v), Wv G 14, 
and the stability estimate for backward Euler schemes gives 

n 

\\8 t u n \\<\\d t u 1 \\ + Y,r\\8 t f J \\. 

1=2 


From (3.3) we have, since Uq = 0, || < ||/i||. Finally, using the inequality 


5^ r ll 5 ‘/lll m ax r ll/(OII <C\\f\\ Lao (L 2 ), 

j =2 1 -^~ 3 


we deduce (3.5). 
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For the bound (3.6) we refer to [HI Theorem 7.3]. 


□ 


3.1. Orthogonal decomposition. In this section we describe the orthogonal de¬ 
composition which defines the GFEM space denoted V in the discussion above. We 
refer to m and [IS] for details. 

For the construction of the GFEM space we use the (weighted) Clement inter¬ 
polation operator introduced in 3J, Jh ■ Vh —> Vh defined by 


3hv = '^2(3 H v)(x)ip x , where (3 H v)(x) := 


Jn v V* 


(3-7) _. . ... 

For this interpolation operator the following result is proved [2] 

(3.8) H^\\v - 3 h v\\l 2{ k) + IIV(v - 3 H v)\\ L2(K) < C\\Vv\\ l ^ k) ,Vv £ V, 

where ook ■= U{K £ Th ■ K D K 0} and C depends on the shape regularity 7 . 

Let V f = {z» £ Vh : 3hv = 0} be the kernel of the Clement interpolation operator 
(3.7). This space contains all fine scale features not resolved by Vh- The space 14 
can then be decomposed into 14 = Vh ® V i , where v £ Vh can be written as a sum 
v = vh + v { , with vh £ Vh, v { £ V 1 , and (vh, v { ) = 0. 

Now define the orthogonal projection R f : Vh —t U f by 

a(R { v,w ) = a(v,w) Vw £ V s , v £ 14. 

Using this projection we define the GFEM space, also referred to as the multiscale 
space, 

U ms := V H - R i V H , 

which leads to another orthogonal decomposition V), = U ms ® U f . Hence any 
function v £ Vh has a unique decomposition v = v ms + v { , with v ms £ V ms and 
v f £ V { , with a(u ms ,u f ) = 0. 

To define a basis for U ms we need to find the projection of the nodal basis 
function tp x £Vh- Let this projection be denoted (j> x , so that cf x £ V { satisfies the 
(global) corrector problem 

(3.9) a(4> x ,w) = a((p x ,w), \/w£V l . 

A basis for the multiscale space U ms is thus given by 

Wx -<t> x : x £ JV}. 

We also introduce the projection R ms : Vh —> V ms , defined by 

(3.10) a(R ms v,w ) = a(v,w), Vw £ U ms , v £ Vh- 

Note that R ms = I — R i . For R ms we have the following lemma, based on the 
results in 0 . 

Lemma 3.2. For the projection R ms in (|3. 10 1 and v £ Vh we have the error bound 


(3.11) 


\v-R ms v\\<CH 2 \\A h v\\, v£V h , 


where C depends on a and 7 , but not on the variations of A. 

Proof. Define the following elliptic auxiliary problem: find z £ Vh such that 
a(z,w) = (v — R ms v,w), Vw £ Vh- 
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In |141 Lemma 3.1] it was proven that the solution to an elliptic equation of the 
form 

a(u,w) = (g,w), Vw £ V h , 

satisfies the error estimate 


|||u - R ms u\\\ < CH\\g\\, 

where C depends on 7 and a, but not on the variations of A. Hence, we have the 
following bound for z , 

p-i? ms ^]|| < CH\\v-R ms v\\. 


Furthermore, we note that v — R ms v £ Vh and 

||t> - R ms i;|| 2 = (t> - R ms v, v - R ms v) = a(z, v - R ms v) 

= a{z - R ms z, v - R ms v) < \\\z - R ms 2 ||| |||v - R ms v |||. 


Now, since a(v, w) = ( AhV , w), we get |||i; — -R ms i;||| < C'U||„4f l u|| and (3.11) follows. 

□ 


In particular, if U n is the solution to (3.3), then (3.11) gives 

| Phfn-B t U n 


I U n - R ms U n II < CH 2 ' 


n > 1 , 


\\d t U n - R ms B t U n || < CH 2 \\P h B t f n - BtdtUn 


n > 2 . 


The result in Lemma 3.2 should be compared with the error of the classical Ritz 
projection Rh : V —> Vh defined by a(RhV,w) = a(v,w ), Vic £ Vh . Using elliptic 
regularity estimates, one achieves 

|| R h v - u|| < Ch 2 \\D 2 v\\ < Ch 2 \\Av\\, 


which is similar to the result in Lemma 3.2 However, in this case, C depends on 
the variations of A , as we noted in the discussion in the beginning of this section. 
This is avoided using the I? ms -projection, since the constant in Lemma [3T2| does not 
depend on the variations of A. 

Now let P ms denote the .^-projection onto U ms and define the corresponding 

eb” 


GFEM to problem (3.3); find U, 
(3.12) 


such that U 0 ms = P ms U 0 and 
(d t U™ s ,v) + a(U™,v) = (/„,«), Vc g U ms , 
for n = 1, ...,7V. Furthermore, we define the operator _4 ms : U ms —> U ms by 
(3.13) (A ms v, w) = a{v, w), Vv, w £ U ms . 


3.2. Localization. Since the corrector problems (3.9) are posed in the fine scale 
space V 1 they are computationally expensive to solve. Moreover, the correctors (j) x 
generally have global support, which destroys the sparsity of the resulting linear 
system (3.12). However, as shown in un, (j) x decays exponentially fast away from 
x. This observation motivates a localization of the corrector problems to smaller 
patches of coarse elements. Here we use a variant presented in [7], which reduces 
the required size of the patches. 

We first define the notion of patches and their sizes. For all I\ £ Th we define 
uJk(K) to be the patch of size k, where 


oj 0 (K) := K, 

u k {K) := U {K £Th ■■ K Ci w k -x{K) ± 0}, k = 1, 2,... 
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Moreover, we define V i {u>k{K)) := {w £ V 1 : supp(w) C Uk(K)}. 

Now define the operator R { K : 14 —> V f by 

/ AV7R l K v ■ Vw = / AS7v ■ V7w, Vv £ Vh, w £ V { , 

Jn Jk 

and note that i? f := J2keT h ^k- We now localize the operator R f K by defining 
R l K ,k : V h -t V i {uj k (K)) through 

[ A\/R l K k v ■ Vw = [ AVv ■ Vw, Vv £ V h , w £ V f (w k {K)), 

Ju k (K) ’ Jk 

and we define R k := Y^keT h ^k k- Hence we can, for each nonnegative integer k. 
define a localized multiscale space 

V k ms ■■= V H - R{V h . 

Here the basis is given by {ip x — <pk,x '■ x £ AT}, where <f> k , x = R^Px is the localized 
version of (f x . The procedure of decomposing 14 into the orthogonal spaces H ms 
and V 1 together with the localization of H ms to H fc m8 is referred to as local orthogonal 
decomposition. 

The following lemma follows from Lemma 3.6 in [7jJ. 

Lemma 3.3. There exists a constant 0 < p < 1 that depends on the contrast (3/a 
such that 

lll^-^||| <Ck d / 2 gL k \\\v\l Vv£V h , 
where C depends on (3, a, and 7 , but not on the variations of A. 

Now let R™ s : 14 —> H™ 8 be the orthogonal projection defined by 

(3.14) a(R™ s v,w) = a(v,w), Vw £ H™ 8 . 

Next lemma is a consequence of Theorem 3.7 in [Tj and estimates the error due to 
the localization procedure. 

Lemma 3.4. For the projection l?™ 8 in ( 3.14| we have the bound 

(3.15) ||v - RTv\\ < C(H + k d > 2 ^fWAhvl Vv e V h . 

Here C depends on [3, a, and 7 , but not on the variations of A. 


Proof. The proof is similar to the proof of Lemma |3.2[ Let z £ Vh be the solution 
to the elliptic dual problem 

a(z,w) = (v — Rjf s v,w), Vw £ Vh, 


which gives 


\v - R^f = (V RTv v R? s v) = a(z - RTz, v - H%*v) 


< IIU - Rl 




It follows from Theorem 3.7 in j7] that there exists a constant C depending on ( 3, a, 
and 7 , such that || 0 -ii^ s z||| < C(H + k d/2 p k )\\v - R™ s v\\, with p as in Lemma 


3.3 


Since ( AhV,w ) = a(v,w) we get |||t> — -R£ 18 d||| < CH\\Ahv\\ and (3.151 follows. □ 
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We are now ready to formulate the localized version of (3.12) by replacing F ms 
by V™. The localized GFEM formulation reads; find Z7“® £ F fe ms such that U™ = 
P^ S U 0 and 


(3.16) 


(d t U£ s n , v) + a(U£ s n , v) = (/„, v), Vv £ F fc ms , 


for n = 1, where P™ s is the L 2 -projection onto V" fc ms . We also define the 


operator : V h 
(3.17) 


ms 

k 


yms 

v k ’ 


a localized version of (3.13), by 


(A” s v, w) = a(v, w), Vv, w £ 14 ms . 

Next lemma states some important properties of the operators Ah, M ms , and 

Ar- 

Lemma 3.5. The operators Ah, -4 ms , and _4“ s are self-adjoint and positive defi¬ 
nite. Furthermore, the following bound hold 

\\(AF)- 1/2 Prf\\<C\\f\\H-i, Vf£L 2 , 

where C depends on a. 

Proof. The fact that the operators are self-adjoint and positive definite follows from 
the assumptions on A (Al). A proof of the bound can be found in [T3]. □ 

We also define the solution operator E™® = ((/ + rA™ s ) _1 ) ra , such that the 
solution to (3.16), with / = 0, can be expressed as UJj'f = Ejf s n U™ q. For this 
operator we have estimates similar to (3.6). 

Lemma 3.6. For l = 0,1, and v £ L 2 , we have 
+—1 


\\d l t E™P? s v\\ < Ct- 1 ||u||, n > l, |||Pn p r«lll < 


< CC 1/2 I 


n> 1, 


where C depends on the constant Ck = sup s> 0 s fe e s , /3, and a. 

Proof. The operator Aff s is self-adjoint, positive definite, and defined on the finite 
dimensional space F fc ms . Thus, there exist a finite number of positive eigenvalues 
{A}fci and corresponding orthogonal eigenvectors {ipi}fL 1 such that span-j^} = 
Vf n!i . We refer to [15] for a further discussion on the eigenvalues to the operator 
A™. 

It follows that E™ s n v can be written as 

M 

E r» = £( 

i—1 x 7 

and the estimates now follows from Lemma 7.3]. 

4. Error analysis 


□ 


In this section we derive error estimates for the local orthogonal decomposition 
method introduced in Section [3| The localized GFEM solution (3.16) is compared 
to the classical FEM solution (3.3), which leads to a setting where the initial data is 
not smooth, since Uq £ 14 only. This leads to error bounds which are non-uniform 
in time, but of optimal order for a fix time t n > 0. The same phenomenon appears 
in classical finite element analysis for equations with nonsmooth initial data, see 
|18j and references therein. The error analysis in this section is carried out by 
only taking the L 2 " n orm of C/ 0 , which allows Uq £ L 2 . If we, for instance, choose 
U 0 = P h u 0 , then ||t/ 0 || < ||w 0 ||. 
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Theorem 4.1. Let U n be the solution to (3.31 and the solution to (3.16|. 
Then, for 1 < n < N, 

lltc -U n \\ < c(l + log ^(H + k^^fit-'WUoW + \\f\\ Loo{L2 ) + ||/IU oo(£a )). 

where C depends on /3, a, 7 , and T, but not on the variations of A. 


The proof of Theorem |4.1| is divided into two lemmas. The first covers the homo¬ 
geneous case, / = 0 , and the second covers the nonhomogeneous case with vanishing 
intital data uq = 0. To study the error in the homogeneous case we use techniques 
similar to the classical finite element analysis of problems with nonsmooth initial 
data, see [15] and the references therein. 

Define XK = ApPh and T™ s = (-4™ S ) _ 1 P™ S . With this notation the solution 


to the parabolic problem (3.3), with / = 0, can be expressed as Thd t U n + U n = 0. 


Similarly, the solution to (3.16), with / = 0, can be expressed as Tff s dtUff*f + 
[/™ = 0. Note that Tjf s is self-adjoint and positive semidefinite on L 2 , and that 

j>ms _ R ms Th 

Now, let e n = t/™® — U n , where e n solves the error equation 

(4.1) T^B t e n +e n = -U n - T? s d t U n = (T h - T? s )5 t U n = (f?^ s - I)U n =: p n , 

for n = 1, ...,N with T™ s e 0 = 0, since /7 ™q = Pfc^Po- The following lemma is a 
discrete versions of m Lemma 3.3]. 

Lemma 4.2. Suppose e n satisfies the error equation (14. It) . Then 


(4.2) \\e n \\ 2 <C(\\p n \\ 2 +t 


n 

i 1 ( 5 Z r ii ft 

1=1 


XX^lH 9 Pi 

i=2 


n > 2 , 


(4-3) INI < || Pl ||. 


Proof. Multiply the error equation (4.1) by d t e n and integrate over S2 to get 


(Xfc St Cn - If e, } ) -f- (e r! . Of e n ) — ( Pm St e n), 

where the first term on the left hand side is nonnegative, since T™ s is positive 
semidefinite on L 2 . Multiplying by rt n we have 


^n||cn|| t n (e n ,e n —\f SL t n (p n ,e n c n _i), 


which gives 


|||e„H 2 - ^||e n - 1 || 2 < t n (p n ,e n - e^) + ^^||e n _ 1 || 2 


— t"n[Pm€"n) ^n— 1 (Pn— 1 5 &n— l) 

(j'nPn tn—lPn—1 5 ^n— 1) H - ^ ||^n—11| • 


Summing over n now gives 


^n||^n|| ^11|^11| — 1 ) ^ ^ J'jPj tj—lPj — li&j—l > ) 

3 = 2 
n 

+ XI r ll e l- 1 l| 2 ’ 

1=2 
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and thus, 

n n 

tn \ \^n || 2 < C(t n \\p n \\ 2 + ^ T {tj\\dtPj\\ 2 + ||Pj —11| 2 ) + ^2 r ll e i-l|| 2 ) • 

3 = 2 3=2 

To estimate the last sum we note that, since T™ s is self-adjoint and positive semi- 
definte, 


2(T’ 


rdtcn 


e n) = (T™ s d t e n , e n ) + {T™ s e n 


i Ui ) 


= dt{T™ s e n , e n ) + T(T£ s 8 t e n , 8 t e n ) > d t (T^ s e n , e n 


so by multiplying the error equation (4.1) by 2e n we get 

<9 t (T“ s e„,e n ) + 2||e n || 2 < 2 (T™ s d t e n ,e n ) + 2||e„|| 2 = 2(p„,e„). 


Multiplying by r and summing over n gives 

n n 

T(T^ s e„,e n ) + ^T|| ej -|| 2 < ^r||pj|| 2 , 
j=i l=i 


where we have used that T™ s eo = 0. Since the first term is nonnegative we deduce 


that X^=i r ll e ll| 2 — Ei=i T ||p/ 1| 2 and (4.2) follows. For n = 1 

(PI). 


this also proves 
□ 


Next Lemma is a discrete version of a result that can be found in the proof of 
[T5T Theorem 3.3]. 

Lemma 4.3. Under the assumptions of Lemma \f.2\ we have, for n > 2, the bound 

j 

(4.4) ||e n || < Ct-\^J\\8 tPj \\ + max^ (ij||pj|| + ||^ r ^ll))- 

— — r—1 

Proof. It follows from Lemma pi~2] that 

||e„|| < C{ max ^||0 t pj|| + max ||py||), n > 2, 

^-SjS n 

or by using Young’s inequality with different constants the proof can be modified 
to show that 


|Gn|| < e max £d|<9 t pd| +C(e) max ||p,||, n > 2, 

2<j<n l<j<n 


for some e > 0. Now define z :] = tje.j . Then 


Tff s d t z n + z n = t n p n + T“ s e n _i := r) n , n > 1, 


and, since = 0 we conclude from Lemma 4.2 


||zj < e max Lj|9 t m|| + C max 

2<j<n 1 <j<n 

From the definition of rjj it follows that 

Ill'll < tj\\pj II + II^T'ej-ill, 3 > 1- 
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Furthermore, for j > 2 

A? II e j —ill — ^'ll^iPjll +tj\\Pj-i\\ i e j— 11 

— ^ll^tPjll +2tj||pj p j—1|| +2tj||pj|| + tj||ej_i|| + 

— 3tj || +2tj||pj|| + 2tj_i 11 e j— 11| 


< 


c(t 2 j\\dtPj\\+tj\\ Pj 


+ 2||zj-i||, 


where we used \tj < tj_ i < tj for j > 2. To bound ||T™ s e n || we define e„ = 
J2j = i Te j and Co = 0. Multiplying the error equation (4.1 ) by t and summing over 
n gives 

n 

TT^ s d t e n +e n = Tjf s d t e n + e n = p n , n > 1, 
l=i 

where p n = Y^j=i T Pj and we have used that T™ s eo = 0. Note that by definition 
T™ s eo = 0. Thus, by Lemma 4.2 we have 


||e„|| < C( max t 3 \\d t pj\\ + max < c( 

\2<j <n / \ 

Hence, since T™ s d t e n = TJf s e n , 


max tj ||ft || + max || Vrp r 

2<j <n 1 <j <n z —' 

r—1 


iTTenll < IISn| 


J 

+ ||Pn|| < C[ max tjUpjll + max ||VTp r ||). 

\ 2 < 7 <n l< 7 <n z —* / 


With e = 4 we get 


\z„\\ < - max + C max H^l 

4 2<j<n l<j<n 


i / - J 

O ma? Ill'll + <4 max ^||0 t pj|| + max (tj||Pjll + II Y' T p, 

/ i < n < m \ 7<n 1 < -»< r> » -* 


< 2 T<i< 


l<;<n 


but from (4.3) we deduce ||zi|| < ti||pi||, and hence 


z n < 


l , i 

a “IK ll^ll + Cl max t|||<9 t ft|| + max (tj||pj|| + || V rp, 

£ 2 7 \^ \ ?< n< n. 1 < -»< r> » ■* 


v 2 <j<n 


l<j<n 


Choosing n* such that ma x 2 <j< n Zj = z n * we conclude (4.4). 


□ 


Lemma 4.4. Assume f = 0 and let £7™ &e the solution to (3.16) and U n the 


solution to (3.3). Then, for 1 < n < N, 

\\uc - u n \\ < c{h + k^pyywuow 

where C depends on /3, a, 7, and T, but not on the variations of A. 


Proof. From Lemma 4.3 we have 


j 

n|| < Ct- 1 (m^J\\B t p j \\ + (tj\\pj\\ + || X] T Pr||)), n > 2, 


r—1 
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and from Lemma 4.2 ||ei|| < ||pi||. The rest of the proof is based on estimates for 
the projection i?™ s in Lemma 3.4 and the regularity of the homogeneous equation 
(|3.6). We have 


t 2 \\d tPj \\ < C(H + k d / 2 n k ) 2 t?\\A h 5 t Uj\\ 

< C{H + k d ' 2 n k ) 2 1? j \\& t 5 t U j \\ < C(H + k d ' 2 p, k ) 2 \\U 0 \\, j > 2 , 
tjWPiWW < C(H + k d ' 2 f i k )H j \\A h U j \\ < C(H + k d / 2 p k ) 2 \\U 0 \\, j > 1, 


E T Pr ii = ii E - T r mil < i w h - irm - u 0 ) 

r= 1 r= 1 

< C(H + k d l 2 n k ) 2 \\Uo\\, 


where we have used ||£/ n || < 11 C/q 11 , which completes the proof. 


□ 


The next lemma concerns the convergence of the inhomogeneous parabolic prob¬ 
lem ( 2.11 with initial data Uq = 0 . 


Lemma 4.5. Assume Uq = 0 and let U™ 8 be the solution to (3.16) and U n the 
solution to (3.3). Then, for 1 <n<N, 


||£4 m : - C7r.ll < C( 1 + log E(iC + fc d /V) 2 (ll/IU oo( L 2 ) 




where C depends on ft, a, 7 , and T, but not on the variations of A. 


Proof. Let U^ n - U n = t / fc m8 - Rff s U n + R^ s U n - U n =: 6 n + p n . For p n we use 
Lemma ITTTtI to achieve the estimate 

\\Pn\\ < C(H + k d / 2 p k ) 2 \\A h U n \\. 

Now, for v £ U^f 118 we have 

(d t 0 n ,v) +a(9 n ,v) = (- Btp n ,v ). 

Using Duhamel’s principle we have 

n 

e n = r E B M-i+i- p r(-3tft), 
i=i 

since 0q = 0. Summation by parts now gives 

n 

e n = E&PTpo - PrPn + r E d t E^ n _ j+1 Pr Pj . 

3 = 1 

Note that po = 0. Using Lemma [3~i] and Lemma [T6| we get 

n 

\K\\ < llPrrll +rJ2 t n-j+l\\Pj\\ 

i=i 

n 

< C(H + k d/2 p k ) 2 max \\A h Uj\\(l + rJ2 t n-j+i)^ 

l<j<n z ' J 

3 =1 

where the last sum can be bounded by 

" t 

t E^-j+i < 1 + logf ■ 

i=i 
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It remains to bound ||„4jiI7 n ||. We have AhU n = Phfn — d t U n and Lemma 
UhUjW < ||/,|| + \\8tUjW < c(||/|| Loo(i2) + \\f\\ Loo(L2) ), 
which completes the proof. 


3.1 


gives 


□ 


Proof of Theorem f.l The result follows from Lemma |4 . 4| and Lemma [4.5| by rewrit¬ 
ing U n = U n< i + U n ^ 2, where [/„,i is the solution to the homogeneous problem and 
U n _2 the solution to the inhomogeneous problem with vanishing initial data. □ 

Remark 4.6. We note that the choice of k and the size of p, determine the rate of 
the convergence. In general, to achieve optimal order convergence rate, k should 
be chosen proportional to log(7L _1 ), i.e. k = c log(ff -1 ). With this choice of k we 
have \\U^ -U n \\< C( 1 + logn)# 2 *” 1 . 

5. The semilinear parabolic equation 

In this section we discuss how the above techniques can be extended to a semi- 
linear parabolic problem with multiscale diffusion coefficient. 

5.1. Problem formulation. We are interested in equations of the form 
u- V- (AVtt) =/(«), in fl x (0, T], 
u = 0, on <9f2 x (0,T], 
u(-, 0) = u o, in 12, 

. is twice continuously differentiable and 12 is a polygonal/polyhedral 
for d < 3. For d = 2,3, / is assumed to fulfill the growth condition 

|/ (,) ( 0 |< C '(1 + |£| 4+1 - | ’> 


(5.1) 


where / : M - 
boundary in 


(5.2) 


), for 2 = 1,2, 


where 6 = 2 if d = 3 and 5 £ [l,oo) if d = 2. Furthermore, we assume that the 
diffusion A fulfills assumption (Al) and uq £ V. 

Example 5.1. The Allen-Cahn equation u — V • (AVu) = —(zi 3 — u) fulfills the 
assumption (5.2). 

Define the ball Br := {c £ V : Mm < R}. Using Holder and Sobolev inequal¬ 
ities the following lemma can be proved, see [T3]. 


Lemma 5.2. If f fulfills assumption ( |5.2| ) and u,v £ Br, then 
\\f(n)\\<C, \\f'(u)z\\ H -i<C\\z\\, \\f'(u)z\\<C\\z\\ H i, \\f(u)z\\ H ^ < C\\z\\, 
and 

\\f(u) - f(v)\\ H -i < C\\u-v\\, 
where C is a constant depending on R. 


From (5.1) we derive the variational form; find u(t) £ V such that 
(5.3) (u, v) + (AVzt, Vi>) = (f(u),v), Vv£V, 

and zt(0) = Uq. For this problem local existence of a solution can be derived given 

that the initial data zto £ V, see 1131 - 


Theorem 5.3. Assume that (Al) and (5.2) holds. Then, for uq £ Br, there exist 
T 0 = T 0 (R) and c > 0, such that (5.3) has a unique solution u £ C(0,T 0 ;U) and 

II m I|l oo (0,T 0 ;U) < cR. 
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For the Allcn-Cahn equation it is possible to find an a priori global bound of 
u. This means that for any time T there exists R such that if u is a solution then 
IKOlUc* ( h 1 ) < R f° r t G [0, X 1 ]. Thus we can apply the local existence theorem 
repeatedly to attain global existence, see m- 


5.2. Numerical approximation. The assumptions and definitions of the families 
of triangulations {77i}/i>o and {Th}h>h and the corresponding spaces Vjj and Vh 
remain the same as in Section [3] For the discretization in time we use a uniform 
time discretization given by 

(5.4) 0 = to < ti < ... < tm = To, where t n — t„_i = r, 

where To is given from Theorem |5.3[ With these discrete spaces we consider the 
semi-implicit backward Euler scheme where U n £ Vf, satisfies 

(5.5) (B t U n , v ) + (AVU n , Vv) = (/(Tn—i), v), to £ V h , 

for n = 1, N where To £ Vi is an approximation of uq. It is proven in [12] that 
this scheme satisfies the bound 


\\U n ~u(t n )\\ <Ct-V\h 2 +T), 

if we choose, for instance, U 0 = PhU 0: where Ph denotes the L 2 -P r ojection onto V) ( . 
Note that C in this bound depends on the variations of A. 

The following theorem gives some regularity estimates of the solution to (5.5). 


Theorem 5.4. Assume that (Al) and (5.2) holds. Then, for Uq £ Bp>, there 
exist T 0 = To(I?) and c > 0 such that (5.5) has a unique solution U n £ Vh, for 
1 < n < N, and maxi<„<Ar ||T n ||jji < cR. Moreover, the following bounds hold 

\\d t U n \\ < Ct~ 1/2 , n > 1, \\\8tU n HI < Ct~\ n > 1, \\B t d t U n \\ < Ct~ 3/2 , n> 2, 

where C depends on a, T 0 , and R, but not on the variations of A. 

Proof. We only prove the estimate \\BtB t U n \\ < Ctf 3 ^ 2 here. The other two follow 
by similar arguments. 

From (5.5) we get 

(5.6) (B t B t U n ,v) + a(B t U n ,v) = (B t f{U n -i),v), Wv£V h ,n> 2 , 

(5.7) {B < f ) U n ,v)+a{B t B t U n ,v) = {B t B t f{U n . l ),v), Wv£V h ,n> 3. 

Choosing v = B t B t U n in ( |5.7| gives 

1 \\B t B t u n \\ 2 - 1 (B t B t u n _ u B t B t u n ) + \\\B t B t u n \f = (B t B t f(u n ^),B t B t u n ), 

T T 


which gives the bound 

( 5 . 8 ) \\B t B t u n \\ 2 - \\B t B t u n -i\\ 2 < CT\\BtBtf(U n -i)\\ H -i. 
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Using Lemma 5.2 we have for £ (min{C/„_j, max{[/„_j, 


||5 t 9 t /( t /„)|| ff - 1 = -||/'(^)( C /„ - U n —\ ) - /'(&)(£/„_! - ^,- 2 )||h-i 

T Z 

< ■4ll(/ , (6)-/ , (6))(^n-Un- 1 )||ir-i 


+ ^2 ll/ / fe)(U n — 2U n -\ + C/n—2) || TT — 1 
< ^II(Ci - C 2 )(t/„ - £/„_!)!! + c||(9 t a t t/ n ||, 

Note that |Ci — £ 2 || < ||17 n _ 2 — £/ n _i|| + ||£/„_i — f/„||. By using Sobolev embeddings 
we get 

4ll(Ci -&){U n - U„_i)|| < max 2||(a t t/,-) 2 || < max 2||0 t [/-||| 

T n—1< j <n n—1<^ <n 

< C max \\d t Uj\\ 2 Hl < Ct\_ x < Ct 2 n , 


where we recall the bounds ^tj < tj _\ < tj for j > 2. Multiplying by rt ^ in (5.8) 
and summing over n gives 

n 

4|| W„|| 2 < ^||W 2 || 2 + J2( Tt 'j\\d t dtf(Un-i)\\ 2 H-i + {t* - tUmdtUj-if) 

3 =3 

n 

< 4\\B t d t u 2 1 | 2 + c^T^imuj^w 2 + +t*_ 1 \\d t dtu j -i\\ 2 ) 

j = 3 

n 

<^p t a t c/ 2 || 2 + ut ri + c^T(^ 4 _ 1 p t 5 t t/ J _ 1 || 2 + t J 3 _ 1 ||a t a t c/,_ 1 || 2 ), 

3 = 3 

for n > 3. Using ||9 t f/j|| < CtJ 1 ^ 2 for j > 1 we get 

4\mu 2 f < CT 2 (||5 t U 2 || 2 + ll^Uill 2 ) < Cr 2 ^- 1 +t- 1 ) < Ct. 


Now, to bound J2j=2^j we choose v = d t d t U n in (5.6) to derive 


(5.9) 


||<W n || 2 + -\\\B t U n III 2 - — \\\3tU n -i ||| 2 < ||5 t /(t/„_ 1 )|| 2 . 


and with C as above, we get 


\\d t f(u n -i)\\ = Wf'i&BtUn-xW < c'ina t ^ n _ 1 ||| < ct-\, 


5.2 


where we used Lemma 

and summing over n gives 


and |||9 t f7j||| < Ct- 1 for j > 1. Multiplying (5.9) with 


^rtf||d t <9 t U^ ; +4|||<9 t Uj| 2 < + (*? “ ^-1)111^^-11 

3 = 2 3 = 2 

+ ^p t c/ 1 |f 


< + rt 2 - i p t C/,- 1 ||| 2 ) +t?p t u 1 ||| 2 . 
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Using |||e? t {7j||| < Ctj 1 for j > 1 we get 

n 

Y r^wBtdtUjf < C(t 2 n + t n + h) < ct n , 

3=2 

where C now depends on t n < T. So we have proved 

n 

t*\\d t 8 t U n \\ 2 < CY rt/ j-i\\dtdtUj-i\\ 2 + Ct n + r 

3=3 

n— 1 n— 1 

^ C 'E rt 3 4 -ill^ f7 t-il| 2 + C'Wi < C'^rt|_ 1 p t a t t/,_ 1 || 2 + Ct rl . 

3=2 3=2 

Applying the classical discrete Gronwall’s lemma gives 

tiWdtdtUjf < Ct n , 

which proves \\BtdtUj\\ < Ctf 3 ^ 2 for n> 3. For n = 2 we proved 

4ll^ t C/ 2 || 2 < CV < Cta, 

which completes the proof. □ 

We use the same GFEM space as in Section [3J that is, U ms = Vh — R { (Vh) and 
the localized version V) ms = Vh — R^Vg)- Furthermore, for the completely dis¬ 
crete scheme, we consider the time discretization defined in (5.41 and the linearized 
backward Euler method thus reads; find U € V ms such that U ™g = Pff^Uo and 

(5.10) m^,v) + a(t/ fe m » = 

for n = 1, ...,7V where P™ s is the .^-projection onto Vf ns ■ 

To derive an error estimates we represent the solution to (5.101 by using Duhamel’s 
principle. Note that E/™® is the solution to the equation 

dtU^ n + ArU^n = PJTf( u k%-i)> 

and by Duhamel’s principle we get 

n 

tt ms 77ms t t ms . \ ^ rms nms unms \ 

0 / fc,n ^k,nT k ,0 ' ^ JUk,j — l)‘ 

3 = 1 

5.3. Error analysis. For the error analysis we need the following generalized dis¬ 
crete Gronwall lemma, see, e.g., Ha- 

Lemma 5.5. Let A, B > 0, 71,72 > 0, 0 < to < t n < T, and 0 <</?„< R. If 

n—1 

Tn < At“ 1+71 +BtJ 2 Y- + i+lTj, 

3=1 

then there is a constant C depending on B, 71 , 72 , and, T, such that, 

in < 43 “ 

Tn V sil n 


Theorem 5.6. For given R > 0 and To > 0 let U n be the solution to (5.5) and 
[7™ be the solution to (5.10), such that U n , £ Br. Then, for 1 < n < TV, 

(5.11) II Uft - U n || < C((H + fc d /y ) 2 + r)t~ 1/2 , 

where C depends on /3, a, 7, R, and To, but not on the variations of A. 
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Proof. First we define e n = U™ s n -U n = - R™ s U n ) + (R™ s U n - U n ) = 9 n +p n 


For pj we use Lemma 3.4 to prove the bounds 

\\ Pj \\<C(H + k d ' 2 p, k )Hj 1/2 , j> 1, 

and 

WdtPkJ < C(H + k d / 2 p k )H~ 3/2 , j > 2. 


For 6 n we have 

n 

On = Ej?J o + r^£^_ j+1 Pr(/(^_i) - / ( Uj — i) - B tPj ). 

7=1 

To bound ||0fe,ra|| we first assume n > 2 and use summation by parts for the first 
part of the sum. Defining 0.2 to be the integer part of nj 2 we can write 


~rJ2Ern-3 + i P rdt Pj = E^PfSpo - E^_ n Pr P n 2 

3 =1 

n 2 

+ t'Z& E ^ 3+i)F™PP 

3=1 

and 9 n can be rewritten as 


r\ pms pms pms 

“n ~ ^k,n^k e 0 &k.n- 


pms 

k,n—ri2 Jr k 


n 2 

p n2 +^E(^- J+ i) p rft 

3=1 


~r E K^j+iPTdtPj 

j—n 2 +1 
n 

+TE(ir) 1/2 ^-i + iWT 1/2 r(/(^-i) - m-i)), 

j=i 

where we note that P™ s eo = 0. To estimate these terms we need the following 
bounds for 71,72 > 0 


_ j-!+7i J-- 1+72 p< J--1+71+72 j--7l j.-1+72 pi 4.-71+72 

' / Pn-i+l L i — *-'7i>72 l 'ra ) ' / Ura-i+lS’ — '- y 7i,72 t n 


7=1 


3=1 


see [12]. Using Lemma 3.6 we get 


11 On 11 < \\pk,n 2 \\ + CT^tn-j+lWPkjW + CT E II d tPk,j 
3=1 j—n 2 +1 

n 

+ ^E*nE+lll/(^b-l) - /(^-l)IU-U 

7=1 
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and together with Lemma [XT] and Lemma [572] this gives 

n -2 n 

iim <c(iL+fc d /V) 2 (c /2 + T E i »Vi i 7 1/2+T E E /2 ) 

i=1 j=n 2 + 1 

n 

+ CrJ2C- / ] 2 + i\\Uk:U-Uj-i\\ 

i=i 

n 

< C(H + fc^V) V /2 + Cr ^C&tllei-rll- 

i=i 

Now consider (?i. We can rewrite 

9 X = EfiOo + rE™P™(f(U^) - /(C7i) - 5 tPl ) 

= ^P fc ms eo + P^P fc m Vi + rE^PnnU^) f(U i)), 

and using similar arguments as above 

l|0i|| < )V /2 + rt7 1/2 ||e 0 ||, 


Hence, we arrive at the estimate 

n 

||e„|| <Ct~ 1/2 {H+ k d/2 fi) 2 +CTj2 t n- J 2 +i\\ e j-i\\> 71 > !. 


j'=i 


and we can use Lemma 5.5 to conclude (5.11). 


□ 


6 . Numerical Results 


In this section we present some numerical results to verify the predicted error 
estimates presented for the linear problem in Section [4] and the semilinear problem 
in Section[5] In both cases the domain is set to the unit square fl = [0,1] x [0,1] and 
T = 1. The domain 11 is discretized with a uniform triangulation and the interval 
[0,T] is divided into subintervals of equal length. 

The method is tested on two different types of diffusion coefficients A\ and A 2 
defined as 


A\{x) 



A 2 {x) 


(B{x) 0 \ 

1 0 B(x)J 


where B is piecewise constant with respect to a uniform Cartesian grid of size 2 -6 , 
see Figure [T] Note that our choice of B imposes significant multiscale behavior 
on the diffusion coefficient. Here we expect quadratic convergence in space of 
the standard finite element with piecewise linear and continuous polynomials (Pl- 
FEM) when A = A 1, but poor convergence when A = A 2 . For the GFEM we 
expect quadratic convergence in both cases. 

We compute the localized GFEM in (3.16) and (5.10), denoted 1/™®, for 5 dif¬ 
ferent values of the coarse grid width, H = y/2 ■ 2 -2 , \/2 ■ 2 -3 , \[2 ■ 2 -4 , \/2 • 2 -5 , 
and %/2 • 2 -6 . The time step is chosen to r = 0.01 for all problems. The reference 
mesh % is of size h = \[2 ■ 2~ 7 and defines the space 14 on which the localized 
corrector problems 4>k, x are solved. To measure the error, the solution U n in (3.3) 
is computed using Pl-FEM on the finest scale h = y/2 ■ 2~ 7 with r = 0.01. 
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(a) Coefficient B for the linear parabolic (b) Coefficient B for the semilinear para- 
problem. The contrast is /3/a « 10 6 . bolic problem. The contrast is /3/a « 10 3 . 


Figure 1. Coefficients for the two parabolic problems. 


Note that this experiment measures the error \\U n — 17™® |. The total error 
\\u{t n ) — 17“® || is also affected by the difference || u(t n ) — U n \\, which is dominating 
for the smaller values of H. We now present the result in two separate sections. 


6.1. Linear parabolic problem. For the linear parabolic problem (2.11 the right 
hand side is set to f(x,t ) = f, which fulfills the assumptions for the required 
regularity. For simplicity the initial data is set to uo = 1. Moreover, at each cell in 
the Cartesian grid we choose a value from the interval [10 _1 ,10 5 ]. This procedure 
gives B a rapidly varying feature and a high contrast max(!?)/min(i3) « 10 6 , see 
Figure [T] (left). 

For each value of H the localized GFEM, 17™®, and the corresponding Pl-FEM, 
denoted Uh,u, are computed. The patch sizes k are chosen such that k ~ log(H” 1 ), 
that is k = 1,2,2, 3, and 4 for the five simulations. When computing Uh, u the 
stiffness matrix is assembled on the fine scale h and then interpolated to the coarser 
scale. This way we avoid quadrature errors. The convergence results for A\ and A 2 
are presented in Figure[2] where the error at the final time fjv is plotted against the 
degrees of freedom |3V"|. Comparing the plots we can see the predicted quadratic 
convergence for the localized GFEM. However, as expected, the Pl-FEM shows 
poor convergence on the coarse grids when the diffusion coefficient has multiscale 
features. We clearly see the pre-asymptotic effects when H does not resolve the 
fine structure of B. 


6.2. Semilinear parabolic problem. For the semilinear problem we study the 
Allen-Cahn equation, which has right hand side f{u) = — (u 3 — u) that fulfills the 
necessary assumptions. We define the initial data to be uo(x, y) = x(l — x)y( 1 — y ), 
which is zero on dfl. The matrix B constructed as in the linear case but with 
values varying between ICC 3 and 1. Note that the solution to the Allen-Cahn 
equation converges to zero rapidly if the diffusion is too high, thus the smaller 
contrast max(H)/min(13) ss 10 3 in this case, see Figure |T] (right). However, B is 
still rapidly varying. As in the linear case we now compute the localized GFEM 
approximations f/™® and the corresponding Pl-FEM, Uh,u- The patch sizes are 
chosen to k = 1, 2, 2, 3, and 4, for the five simulations. The convergence results for 
Ai and A 2 are presented in Figure [3] We can draw the same conclusions as in the 
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Figure 2. Relative L 2 errors \\U™lf — Uh,N\\/\\Uh,N\\ (blue o) and 
\\Uh,n — U hiN \\/\\U htN \\ (red *) for the linear parabolic problem 
plotted against the number of degrees of freedom \AF\ ~ H~ 2 . The 
dashed line is H 2 . 

linear case. The localized GFEM shows predicted quadratic convergence in both 
cases, but Pl-FEM shows poor convergence on the coarse grids when the diffusion 
coefficient has multiscale features. 



Figure 3. Relative L 2 errors — {7/ 1i jv||/||?7/i,jv|| (blue o) and 
\\UH,N—Uh,N\\/\\Uh } N\\ (red *) for the semilinear parabolic problem 
plotted against the number of degrees of freedom \AF\ « H~ 2 . The 
dashed line is H 2 . 
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